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The invention concerns an image metiiod for observing the propagation of low-frequency shearing pulse 
wave simultaneously in multiple points of a diffusing viscoelastic medium! The method consists in 
transmitting at a very high rate ultrasonic compression waves enabling to obtain a succession of images 
of the medium; then in delayed processing of the resulting images by intercorrelation to determine in 
each point of each image the movements of the medium while the shearing wave is being propagated. 
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PROCEDE ET DiSPOSITIF D'IMAGERIE UTiLISANT LES ONDES DE CiSAILLEIUIENT. 

(57) Procecl6 dNmagerie pour observer la propagation 
d une onde impulsionnelle de cisaillement basse frequence 
simuttan^ment en une multitude de points d'un milieu vis- 
co^lastique diffusant (1). A cet effet, on 6met k cadence ul- 
trarapide des ondes uitrasonores de compression qui 
permettent d'obtenir une succession dMmages du milieu, 
puis on tralte en temps diff6r6 les images ainsi obtenues par 
mtercorr^lation, pour determiner en chaque point de cheque 
Image les mouvements du milieu lors de la propagation de 
I'onde de cisaillement. 
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p-rocede et disyr^flitiif r! ' Imaqrerie utilisant — 1 ^ff On^gg 
cisaillement . 

La presente invention est relative aux precedes et 
dispositifs d' imager ie ut 11 isant les ondes de cisaille- 
ment . 

Plus particuliSrement, 1» invention concerne un 
precede d'imagerie utilisant les ondes de cisaillement 
pour observer un milieu viscoelastique diffusant qui con- 
tient des particules ref Ifichissant les ondes ultrasonores 
de compression, proced§ dans lequel on genere une onde 
elastique de cisaillement dans le milieu viscoelastique et 
on observe, au moyen d'au moins une onde ultrasonore de 
compression, le deplacement du milieu viscoelastique sou- 
mis a ladite onde de cisaillement. 

Le document US-A-5 810 731 decrit un exemple d'un 
tel precede, dans lequel 1 ' onde de cisaillement est genS- 
ree localement a 1 ' interieur du milieu viscoelastique ob- 
serve, au moyen de la pression de radiation d'une onde ul- 
trasonore modulee et focalisee sur un point a observer. On 
envoie ensuite sur ce point focal une onde ultrasonore 
supplementaire dont la reflexion permet de connaitre cer- 
tains parametres de propagation de 1 • onde de cisaillement 
(notamment la viscosite dynamique du milieu ec son module 
de cisaillement) au niveau du point focal susment ionne . 

Cette technique presente 1 » inconvenient de permet- 
tre 1' analyse d'un seul point du milieu viscoelastique 
etudie a chaque fois qu'une onde de cisaillement est gene- 
ree. Si 1 ' on veut obtenir une image complete du milieu 
viscoelastique observg, il est necessaire de repeter 
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1 'operation un tres grand nombre de fois, ce qui implique 
un temps de pause important (par exemple, plusieurs minu- 
tes) pour obtenir cette image. 

Ce temps de pause important rend ce procede de 
I'art anterieur peu pratique a utiliser. 

De plus, un tel temps de pause peut nuire a 1' uti- 
lisation dudit precede pour obtenir une image d'un tissu 
vivant, qui est toujours en mouvement . 

La prSsente invention a notamment pour but de pal- 
lier ces inconvSnients . 

A cet effet, salon !• invention, un procede du 
genre en question est essent iellement caracterise en ce 
qu'on genere l*onde de cisaillement en appliquant au mi- 
lieu viscoSlastique une excitation ayant la forme d'une 
impulsion basse frSquence qui presente une frequence cen- 
trale f comprise entre 20 et 5000 Hz, cette impulsion 
basse frequence prSsentant une duree comprise entre l/2f 
et 20/f, 

en ce qu ' il comporte une etape d' observation de propaga- 
tion au cours de laquelle on observe la propagation de 
I'onde de cisaillement simultanement en une multitude de 
points dans le milieu observe, ces points formant un champ 
d' observation sensiblement continu s ' Stendant au moins se- 
lon un premier axe, cette etape d ' observation de propaga- 
tion de I'onde de cisaillement consistant a : 

emettre dans le milieu observe une succession 
d'au moins 10 tirs d'ondes ultrasonores de compression a 
une cadence comprise entre 100 et 100 000 tirs par se- 
conde , 
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detecter et enregistrer en temps reel les echos 
generes par les particules ref lechissantes du milieu vis- 
coelastique a chaque tir d'onde ultrasonore, ces echos 
correspondant (directement ou indirectement ) a des images 

5 successives du milieu observe, 

et en ce que ledit precede comporte en outre une ecape ul- 
terieure de traitement d' image au cours de laquelle on 
traice en temps differe les images ainsi obtenues au moins 
par intercorrSlation entre images successives, pour deter- 

10 miner en chaque point du champ d' observation un paramStre 
de mouvement choisi entre le deplacement et la deformation 
du milieu viscoelastique , de fagon a obtenir ainsi une 
succession d' images montrant 1 'Evolution du parametre de 
mouvement du milieu viscoSlastique sous I'effet de la pro- 

15 pagation de I'onde de cisaillement . 

Grace a ces dispositions, on obtient un film il- 
lustrant clairement la propagation de 1 ' onde de cisaille- 
ment dans le milieu viscoelast ique , qui peut permettre par 
exemple, dans des applications medicales, de reperer di- 

20 rectement des zones cancSreuses dans les tissus d'un pa- 
tient la propagation des ondes de cisaillement s'y de- 
roule en effet trSs diffSremment des zones voisines. 

Ce reperage s'effectue beaucoup plus facilement 
que par une observation classique par simple echographie 

25 ultrasonore, puisque la propagation des ondes de cisaille- 
ment est fonction du module de cisaillement du milieu, 
lui-meme tres variable entre une zone de tissus sains et 
une zone de tissus cancereux : le module de cisaillement 
varie typiquement dans un rapport de 1 a 30 entre une zone 

30 saine et une zone cancereuse, alors que le module de com- 
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pression, qui regit la propagation des ondes acoustiques 
de compression utilisees dans 1 • echographie ultrasonore, 
varie seulement de 1 ' ordre de 5% entre un tissu sain et un 
tissu canc§reux. 

On notera que le film obtenu illustre la propaga- 
tion de I'onde de cisaillement beaucoup plus clairement 
que la simple succession des images donnees par les parti- 
cules ref l§chissantes du milieu, puisque ledit film permet 
de visual iser chaque instant les zones du milieu observe 
qui subissent des mouvements de meme ampleur du fait de la 
propagation de I'onde de cisaillement / alors que la suc- 
cession des images des particules ref lechissantes ne per- 
nettrait de visualiser qu'un brouillard de points lumineux 
en mouvement . 

Dans des modes de realisation preferes du precede 
selon invention, on peut eventuellement avoir recours en 
outre a I'une et/ou S 1* autre des dispositions suivantes : 

la duree de 1' impulsion basse frequence est 
comprise entre l/2f et 2f ; 

la frequence centrale de 1' impulsion basse fre- 
quence est comprise entre 30 et 1000 Hz ; 

le milieu viscoelast ique observe est constitue 
par un corps vivant comprenant au moins un organe interne 
soumis a des mouvements impulsionnels , 1' impulsion basse 
frequence qui gendre I'onde de cisaillement etant consti- 
tuee par un mouvement impulsionnel dudit organe interne ; 

le milieu viscoelast ique observe est delimite 
par une surface exterieure, et 1' impulsion basse frequence 
est appliquee au niveau de cette surface exterieure ; 



T 1 
1 

2791136 



5 

1' impulsion basse frequence est appliquee par 
un mcyen d • excitation choisi parmi : 

une onde acoustique generee par au moins un trans - 

ducteur acoustique, 
5 . et un choc gSnSrfi localement par contact physique 

au niveau de la surface exterieure du milieu viscoelasti- 
que ; 

on emet les tirs d'ondes ultrasonores de com- 
pression et on detecte les 6chos generes par les particu- 

10 les ref lechissantes du milieu visco§last ique , au moyen 
d'une batterie de transducteurs qui comprend au moins un 
transducteur et qui est disposee au contact de la surface 
exterieure du milieu visco^lastique , I'onde de cisaille- 
mant etant appliquee au milieu viscoSlast ique en imposant 

15 un deplacement impulsionnel a ladite batterie de transduc- 
teurs ; 

ledit parametre de mouvement est la deformation 
du milieu viscoelast ique : cette disposition est particu- 
lierement utile dans le dernier cas envisage ci-dessus, 
20 puisqu'elle permet de s'affranchir du deplacement de la 
batterie de transducteurs, deplacement qui perturberait 
sinon la mesure du deplacement des points du champ d' ob- 
servation ; 

au cours de 1 ' observation de la propagation de 
25 I'onde de cisaillement , on emet entre 100 et 10 000 tirs 
d'ondes ultrasonores de compression a une cadence comprise 
entre 100 et 100 000 tirs par seconde ; 

le champ d • observation s'etend au m.oins selon 
un plan comprenant d'une part, le premier axe et d' autre 
30 part, un deuxieme axe perpendiculaire au premier axe ; 
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au cours de I'etape d* observation de propaga- 
tion, on utilise une batterie de plusieurs transducteurs 
acoustiques disposes au moins selon le deuxieme axe pour 
emettre les tirs d'ondes ultrasonores de compression et 
5 detecter les echos gSnSrSs par les particules refl^chis- 
sances du milieu viscoelastique, les 6chos decectSs par 
chaque transducteur acoustique Stant memorises directement 
sans traitement prSalable au cours de I'etape d'observa- 
tion de propagation, et I'Stape de traitement d' image com- 

10 prenant une sous-etape prSliminaire de formation de voies 
au cours de laquelle on g^nSre une image du milieu visco- 
elascique correspondant Sl chaque tir d'onde ultrasonore de 
compression, par combinaison d*au moins certains des echos 
re<;:us par les differents transducteurs ; 

15 - I'etape de traitement d' images est suivie (im- 

mediatement ou non) par une etape de visualisation au 
cours de laquelle on visualise au ralenti un film consti- 
tue par la succession des images traitees, chaque point de 
chaque image presentant un parametre optique qui varie 

20 suivant la valeur du parametre de mouvement affecte a ce 
point ; 

le paramdtre optique est choisi parmi le niveau 
de gris et le niveau chromatiq[ue ; 

I'Stape de traitement d' images est suivie (im- 
25 mediatement ou non) par une Stape de cartographie au cours 
de laquelle, S partir de 1' evolution du paramdtre de mou- 
vement au cours du temps dans le champ d • observation, on 
calcule au moins un parametre de propagation de I'onde de 
cisaillement en au moins certains points du champ d'obser- 
30 vat ion ; 
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le parametre de propagation de 1 ' onde de ci- 
saillement qui est calcule au cours de I'etape de carto- 
graphie est choisi parmi la vitesse des ondes de cisaille- 
(nent, le module de cisaillement , 1 ' attenuat ion des ondes 

5 de cisaillement , I'elasticite de cisaillement , et la vis- 
cosite de cisaillement. 

Par ailleurs, 1' invention a egalement pour objet 
un dispositif d'imagerie utilisant les ondes de cisaille- 
ment pour observer un milieu viscoelast ique diffusant qui 

10 contient des particules ref lechissant les ondes ultrasono- 
res de compression, ce dispositif comprenant des moyens 
d' excitation pour gSnfirer une onde elastique de cisaille- 
ment dans le milieu viscoelastique et des moyens d' acqui- 
sition pour observer, au moyen dVau moins une onde ultra- 

15 sonore de compression, le deplacement du milieu viscoelas- 
tique soumis a ladite onde de cisaillement, 

caraccerise en ce que les moyens d' excitation sont adaptes 
pour appliquer au milieu viscoelastique une excitation 
ayanc la forme d'une impulsion basse frequence qui pre- 
20 sente une frequence centrale f comprise entre 20 et 5000 
Hz, cecte impulsion basse frequence presentant une duree 
comprise entre l/2f et 20/f, 

en ce que les moyens d ' acquisition sont adaptes pour ob- 
server la propagation de I'onde de cisaillement simultane- 

25 ment en une multitude de points dans le milieu observe, 
ces points formant un champ d ' observation sensiblement 
continu s ' etendant au moins selon un premier axe, lesdits 
moyens d ' acquisition etant adaptes pour : 

emettre dans le milieu observe une succession 

30 d'au moins 10 tirs d' ondes ultrasonores de compression a 
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une cadence comprise entre 100 et 100 000 tirs par se- 
conde , 

detecter et enregistrer en temps reel les echcs 
generes par les particules ref ISchissantes du milieu vts- 
coelastique k chaque tir d'onde ultrasonore, ces echos 
correspondant (directement ou indirectement ) a des images 
successives du milieu observfi, 

et en ce que ledit dispositif comporte en outre des mcyens 
de traitement d* image adaptSs pour traiter en temps dif fe- 
re les images obtenues Sl partir des moyens d ' observation, 
au moins par intercorrelat ion entre images successives, 
pour determiner en chaque point du champ d' observation un 
paramecre de mouvement choisi entre le diplacement et la 
deformation du milieu viscoelastique , de fagon a obtenir 
ainsi une succession d' images montrant 1' evolution du pa- 
rametre de mouvement du milieu viscoelas t ique sous I'effet 
de la propagation de I'onde de cisaillement . 

Dans des modes de realisation preferes du disposi- 
tif selcn 1' invention, on peut eventuellement avoir re- 
cours en outre a 1 • une et/ou a 1' autre des dispositions 
suivantes : 

les moyens d ' observation comprennent une batte- 
rie de transducteurs qui inclut au moins un transducteur 
et qui est adaptee pour etre disposee au contact d'une 
surface exterieure dSlimitant le milieu viscoelastique , 
les moyens d' excitation etant adaptes pour imposer un de- 
placement impulsionnel Sl ladite batterie de transduc- 
teurs ; 

ledit parametre de mouvement est la deformation 
. du milieu viscoelastique ; 
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le champ d* observation s'ecend au moins selon 
un plan comprenant d'une part, le premier axe et d' autre 
part, un deuxieme axe perpendiculaire au premier axe, la 
batterie de transducteurs comportant plusieurs transduc- 
5 teurs disposes au moins selon le deuxieme axe, des moyens 
de commande etant prevus pour select ivement faire fonc- 
tionner le dispositif soit dans le mode d* imagerie par on- 
des de cisaillement , soit dans un mode d • echographie stan- 
dard permettant d'acquerir entre 10 et 100 images par se- 
10 conde . 

D^autres caracteristiques et avantages de 1' inven- 
tion apparaitront au cours de la description suivante de 
plusieurs de ses formes de realisation, donnees a cicre 
d' examples non limitatifs, en regard des dessins joints. 
15 Sur les dessins : 

la figure 1 est une vue schematique d'un dispo- 
sitif d' imagerie par ondes de cisaillement selon une forme 
de realisation de 1* invention, 

et la figure 2 est une vue de detail montrant 
20 une variante du dispositif de la figure 1. 

Sur les differentes figures, les memes references 
designent des elements identiques ou similaires. 

La figure 1 represente un exemple de dispositif 
d' imagerie par ondes de cisaillement selon 1' invention, 
25 pour etudier la propagation des ondes elastiques de ci- 
saillement dans un milieu viscoelastique 1 qui est diffu- 
sant vis ^ vis des ondes ultrasonores de compression, et 
qui peut etre par exemple : 
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un corps inerte, notamment dans le cas du con- 
trole de qualite pour des applications industrielles, no- 
tammenc agro-alimentaires , 

ou un corps vivant, par exemple une parcie du 
corps d'un patient, dans le cas des applications medica- 
les . 

Ce dispositif comporte un transducteur acoustique 
tel qu'un haut-parleur 2 ou un pot vibreur qui est dispose 
contra la surface extSrieure 3 du milieu observe 1, cette 
surface 3 etant constitute par exemple par la peau du pa- 
tient dans les applications medicales. 

Le haut-parleur 2 peut etre commande par un micro- 
ordinateur 4, par exemple par 1 • intermSdiaire d*un circuit 
generateur d' impulsions basse frequence G (ce circuit peut 
acre constitue notamment par la carte son du micro- 
ordinateur 4) et d ' un amplif icateur A, pour appliquer a la 
surface 3 du milieu observe une excitation sous la forme 
d'une impulsion basse frequence, de fagon a generer une 
onde de cisaillement dans le milieu viscoelast ique 1. 

Cette impulsion basse frequence presente en gene- 
ral une amplitude pouvant etre de I'ordre de 1 mm et une 
frequence centrale f comprise entre 20 et 5000 Hz, appli- 
quee pendant une duree comprise entre l/2f et 20/f. De 
preference, la duree d ' application de 1 ' impulsion basse 
frequence est comprise entre l/2f et 2f et la frequence f 
est comprise entre 30 et 1000 Hz, cette frequence pouvant 
typiquement §tre de I'ordre de 50 Hz. 

En variante I'onde acoustique de cisaillement 
pourrait egalement §tre obtenue (avec les caracttrist iques 
susmentionntes d' amplitude et de frequence) : 
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par un choc genere localement par contact phy- 
sique au niveau de la surface exterieure du milieu visco- 
elastique, par au moins un actionneur mecanique automatic 
que commande par le micro-ordinateur 4, 

par un choc applique manuellement par contact 
physique au niveau de la surface exterieure du milieu vis- 
coelast ique, 

ou encore, par un mouvement impulsionnel natu- 
rel d'un organe interne du corps humain ou animal (par 
exeTiple un battement du coeur) dans les applications medi- 
cales . 

L'onde elastique de cisaillement produite par le 
haut-parleur 2 se deplace avec une vitesse Cs relativement 
faible, de 1 ' ordre de quelques m/s { typiquement , de 1 a 10 
m/s dans le corps humain) en produisant des mouvements in- 
ternes dans le milieu viscoelast iqfue 1 observe. 

Ces mouvements sont suivis en envoyant dans le mi- 
lieu 1, des ondes ultrasonores de compression qui inter- 
agissent avec les particules diffusantes 5 contenues dans 
le milieu 1, lesquelles particules sont ref lechissantes 
pour les ondes ultrasonores de compression. Les particules 
5 peuvent etre constituees par toute heterogeneite du mi- 
lieu 1, et notamment, lorsqu'il s'agit d'une application 
medicale, par des particules de collagene presentes dans 
les tissus humains. 

Pour observer la propagation de l'onde de ci- 
saillement, on utilise done une sonde ultrasonore 6 dispo- 
see centre la surface exterieure 3 du milieu observe 1. 
Cette sonde envoie, selon un axe X, des impulsions d* ondes 
ultrasonores de compression du type de celles couramment 
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utilisees en echographie, a une frequence comprise par 
exemple entre 1 et 100 MHz et de preference entre 3 et 15 
MHz. On nocera que la sonde 6 peut etre disposee : 

soit du meme cSte du milieu 1 que le haut- 
5 parleur 2, comme represents sur la figure 1, 

soit k 1* oppose du haut -parleur 2 par rapport 
au milieu 1, 

soit dans toute autre position, par exeir.ple 
dans une disposition transversale par rapport au haut- 
10 parleur 2. 

La sonde ultrasonbre 6 est coiistituee par une bat- 
ter ie de n transducteurs ultrasonores Tl, T2 , Ti, 
Tn, n etant un nombre entier n au moins egal a 1. 

Cette sonde 6 se pr^sente le plus souvent sous la 
15 forme d'une barrette linSaire pouvant comprendre par exem- 
ple n - 128 transducteurs alignes selon un axe Y perpendi- 
culaire a 1 ' axe X, envoyant simultanement leurs impulsions 
d'ondes ultrasonores de fa^on a generer une onde "plane" 
(c*est a dire en 1* occurrence une onde dont le front 
20 d'onde est rectiligne dans le plan X, Y) ou tout autre 
type d'onde eclairant 1' ensemble du champ d' observation . 

En variante, la batterie 2 de transducteurs peut 
eventuellement etre reduite ^ un transducteur unique Tl, 
ou au contraire se presenter sous la forme d'un reseau a 
25 deux dimensions s'etendant par exemple selon un plan per- 
pendiculaire a 1 ' axe X. 

Chacun des transducteurs Tl, T2, ... Tn est com- 
mande par le micro-ordinateur 4 ou par une unite centrale 
CPU (qui est contenue par exemple dans une bale electronic 
30 que 7 reliee par un cSble souple ^ la sonde 6), pour emet- 
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tre des tirs successifs d'ondes de compression ultrascno- 
res dans le milieu 2, au cours d'une phase d ' observation 
qui peut durer par exemple moins d'une seconde et au cours 
de laquelle sont emis p tirs d'ondes ultrasonores de com- 
5 pression (p ecant un entier compris entre 100 et 10 000 et 
de preference entre 1000 et 100 000) , ^ une cadence com- 
prise entre 100 et 100 000 tirs par seconde et de prefe- 
rence comprise entre 1000 et 100 000 tirs par seconde, no- 
tamment entre 1000 et 10 000 tirs par seconde (cette ca- 

10 dence est limitee par le temps d ' aller-retour de I'onde de 
compression dans le milieu 1, done par I'epaisseur du mi- 
lieu 1 dans la direction X : il faut en effet que tous lea 
echos generes par 1 ' onde de compression aient: ece recpus 
par la sonde 6 avant d'envoyer une nouvelle onde de com- 

15 pression) . 

Les tirs d'ondes ultrasonores de compression de la 
phase d ' observation commencent de preference juste avant 
!• emission de 1 ' onde de cisaillement . 

De plus, dans le cas oCt I'onde de cisaillement est 

20 generee par un mouvement impulsionnel d'un organe d'un 
corps vivant, on peut avantageusement synchroniser le de- 
marrage des tirs d'ondes ultrasonores de compression avec 
ce mouvement impulsionnel. par exemple dans le cas ou 
I'onde de cisaillement est generSe par un battement de 

25 ccBur, on. peut synchroniser le demarrage des tirs d'ondes 
ultrasonores de compression avec une phase choisie de 
i ' electrocardiogramme . 

Chacun de ces tirs donne lieu a la propagation 
* d'une onde ultrasonore de compression dans le milieu 1, 

30 avec une vitesse de propagation beaucoup plus elevee que 
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les ondes de cisaillement , par exemple de 1 ' ordre de 1500 
m/s dans le corps humain. 

L'onde ultrasonore ainsi generee interagit avec 
les particules ref lechissantes 5, ce qui genere des echos 

5 ou aucres perturbations analogues du signal, connus en soi 
sous le nom de "bruits de speckle" dans le domaine de 
1 • echographie . 

Ces "bruits de speckle" sont captes par les trans- 
ducteurs Tl, ,Tn aprSs chaque tir. Le signal sij (t) ain- 

10 si capte par chaque transducteur Ti aprSs le tir n*» j est 
tout d'abord echantillonng Sl haute frequence (par exemple 
de 30 a 100 MHz) et numerise en temps reel (par exemple 
sur 8 bits, ou dans certains cas sur 1 bit) par un 6chan- 
tillonneur appartenant Sl la bale 7 et relie a ce transduc- 

15 teur, respect ivement El, E2 , •» En. 

Le signal sij (t) ainsi echantillonn§ et numerise 
est ensuite memorisS, egalement en temps rSel, dans une 
memoire Mi appartenant ^ la bale 7 et propre au transduc- 
teur Ti. 

20 Chaque mSmoire Mi presente par exemple une capaci- 

te de I'ordre de 1 Mo, et contient 1* ensemble des signaux 
sij (t) regus successivement pour les tirs j = 1 a p. 

En temps differe, apres la memorisation de tous 
les signaux sij (t) correspondant a une meme propagation 

25 d'onde de cisaillement, 1' unite centrale CPU fait retrai- 
ter ces signaux par un circuit sommaceur S appartenant a 
la bale 7 (ou bien elle effectue elle-meme ce traitement, 
ou encore ledit traitement peut etre effectue dans le mi- 
cro-ordinateur 4) , par un processus classique de formation 

30 de voies . 
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On genere ainsi des signaux Sj (x, y) qui corres- 
pondent chacun a 1' image du champ d' observation apres le 
tir n" j . 

Par exemple, on peut determiner un signal Sj (t) 
5 par la formule suivante : 

SjV) = Zor. (X, y)MMx. y) + d, (x,>) / V] 

OU : 

sij est le signal brut pergu par le transduc- 
ceur i apres le tir d'onde ultrasonore de compression 

10 n<* j, 

t(x,y) est le temps mis par I'onde ultrasonore 
de compression pour atteindre le point du champ 
d' observation de coordonnees (x,y), avec t =^ 0 au debut du 
tir n<^ j , 

15 - di(x,y) est la distance entre le point du champ 

d' observation de coordonnees (x,y) et le transducteur n<> 
i, ou une approximation de cette distance, 

V est la Vitesse moyenne de propagation des on- 
des acoustiques ultrasonores de compression dans le milieu 
20 viscoelastique observ§. 

et ai{x,y) est un coefficient de pondSration 
tenant compte de lois d' apodisation (en pratique, on pour- 
ra dans de nombreux cas considerer que ai(x,y)=l) . 

La formule ci-dessus s' applique mutatis mutandis 
25 lorsque le champ d' observation est a 3 dimensions (reseau 
plan de t ransducteurs) , en remplagant les coordonnees spa- 
tiales (x,y) par (x,y,z). 
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Lorsque la sonde 6 comporte un seul transducteur , 
I'etape de formation de voies est inutile, et on a direc- 
tement Sj (x) = sj[2.x/V], avec les memes notations que 
ci-dessus . 

5 Aprds l»etape eventual le de formation de voies, 

1* unite centrale CPU memorise dans une memoire centrale M 
appartenant a la bale 7 les signaux d' images Sj(x,y) ou 
Sj (x) ou Sj(x,y,z), qui correspondent chacun au tir n^ j. 
Ces signaux peuvent egalement etre memorises dans le mi- 

10 cro-ordinateur 4 lorsqu'il effectue lui-meme le traitement 
d ' image . 

Ces images sont ensuite trait^es deux Sl deux, ou- 
jours en temps diffSr^, par intercorrelation. Cette inter- 
correlation peut Stre rSalisSe dans un circuit DSP appar- 

15 tenant a la bale 7, ou Stre programmee dans l'unit6 cen- 
trale CPU ou dans le micro-ordinateur 4 . 

A titre d'exemple, cette intercorrelation peut se 
faire en comparant les signaux Sj (x,y) et Sj+1 (x,y) 
(dans le cas d'un champ d' observation a 2 dimensions) sur 

20 des fenetres spatiales glissantes de longueur Ax prSdSter- 
minee, pouvant aller par exemple de A. a lOA., ou A. est la 
longueur d'onde des ondes ultrasonores de compression 
(soit environ 0,42 Sl 4,2 mm a 3,5 MHz dans I'eau ou le 
corps humain) . De plus, les fenetres susment ionnees peu- 

25 vent se recouvrir I'une 1' autre sur environ 20% de leur 
longueur selon I'axe X. 

Au cours de ce processus d' intercorrelation, on 
maximise une fonction d ' intercorrelation 

<Sj (x,y) ,Sj+l (x,y) > afin de determiner le deplacement subi 
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par chaque particule 5 donnant lieu a un echo ultrasonore, 
dans la direction X. 

Des exemples de tels calculs d' intercorrelat ion 
sont donnes dans I'etat de la technique, notamment par 
5 O'Donnell et al . (« Internal displacement and strain ima- 
ging using speckle tracking », IEEE transactions on ultra- 
sonic, f erroelectrics, and frequency control, vol. 41, 
3, mai 1994, p. 314-325) et par Ophir et al . 
{« Elastography : a quantitative method for imaging the 

10 elasticity of biological tissues », Ultrasonic imag., vol. 
13, p. 111-134, 1991) . 

On obtient ainsi une succession de champs de d6- 
placements Djx(x,y) du milieu 1 sous I'effet de I'onde de 
cisaillement , dans la direction X. 

15 Cette succession de champs de deplacements est 

stockSe dans la memoire M ou dans le micro-ordinateur 4 et 
peut etre visualisee, notamment au moyen de I'ecran 4a du 
micro-ordinateur, sous la forme d'un film ralenti ou la 
valeur des deplacements est illustree par un parametre op- 

20 tique tel que par un niveau de gris ou par un niveau chro- 
mat ique . 

On visualise ainsi parfaitement les differences de 
propagation de 1 * onde de cisaillement entre les zones de 
caracteristiques differentes du milieu 1, par exemple les 
25 tissus sains et les tissus cancereux dans le cas d'une ap- 
plication m^dicale. 

Ce film de propagation de I'onde de cisaillement 
est en outre superposable avec une image echographique 
classique, qui peut §tre ggnSree par le dispositif decrit 
30 ci-dessus, lequel est apte a fonctionner : 
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soit dans le mode d'imagerie par ondes de ci- 

saillement , 

soit dans un mode d' echographie standard, en 
fonction des commandes recjues par example depuis le cla- 
5 vier 4b du micro -ordinateur . 

Par ailleurs, il est ggalement possible de calcu- 
ler non pas les dSplacements de chaque point du milieu ob- 
serve 1, mais les deformations Ejx(x,y) du milieu dans la 
direction X, c'est a dire les derivges des deplacements 
10 Djx(x,y) par rapport a x, 

Ces champs de deformations successifs sont utili- 
sables comme precedemment pour visualiser . clairement la 
propagation de 1 ' onde de cisaillement sous la forme d'un 
film, et presentent en outre 1 ' avantage de s ' af f ranchir ■ 
15 des deplacements de la sonde 6 par rapport au milieu ob- 
serve 1 . 

Cette variants est part iculierement interessante 
dans la forme de realisation de la figure 2, ou la sonde 6 
est portSe par le haut-parleur ou vibreur 2, ce qui impli- 
20 que forcement des mouvements de ladite sonde puisque c'est 
elle-mSme qui genSre alors I'onde de cisaillement. 

A partir des champs de deplacements ou de deforma- 
tions, on peut le cas ^cheant procSder en outre a une 
Stape de cartographic au cours de laquelle, k partir de 
25 1' evolution du parametre de mouvement (deplacement ou de- 
formation) au cours du temps dans le champ d* observation 
X, Y (ou X dans la ces d'un seul transducteur , ou X, Y, Z 
dans le cas d'un reseau plan de transducteurs) , on calcule 
» au moins un param&tre de propagation de I'onde de ci- 
30 saillement, soit en certains points du champ d • observation 



BNSOOCID:<FR 27911 36A1 I > 



2791136 



19 

choisis par 1 ' utilisateur S partir du micro- ordinateur 4, 
soit dans tout le champ d ' observation . 

Le parametre de propagation de I'onde de cisaille- 
ment qui est calcule au cours de I'etape de cartographie 

5 est choisi par exemple parmi la vitesse Cs des ondes de 
cisaillement , le module de cisaillement \i, 1 ' attenuation a 
des ondes de cisaillement, 1' elasticity de cisaillement 
lalviscosite de cisaillement ji2 , 

Ce calcul est effectue par un processus classique 

10 d* inversion, dont un exemple est donne ci-apres dans le 
cas d*un champ d ' observation a deux dimensions (le meme 
processus s ' appliquerait mutatis mutandis dans le cas d'un 
champ d' observation a une ou a trois dimensions, respec- 
tivement pour un seul transducteur Tl ou pour un reseau 

15 plan de transducteurs) . 

Dans cet exemple on fera 1 ' approximation cjue la 
viscosite de cisaillement jj.2 est nulle et que le milieu 
est isotrope- 

L' equation d*onde donnant le vecteur deplacement D 
20 de chaque point du milieu 1 s'ecrit : 

P^ = A(//.D) (I), 
ct' 

ou p est la density du milieu 1, ^ est le module de ci- 
saillement (par hypothdse reduit a sa partie reelle 
I'elasticite de cisaillement |nl puisque la viscosite de 
25 cisaillement ^2 est supposee nulle) . 

Pour la premiere composante u du vecteur D, c'est 
a dire pour le deplacement du milieu 1 dans la direction 
X, on a done : 
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^ dr at- dy"- 

AprSs transformSe de Fourier temporelle et discre- 
tisation, cette equation peut etre ecrite sous, la forme 
matricielle suivante, qui peut etre §crite pour chaque 
5 frequence du spectre de I'onde de cisaillement : 

B^HM (III) / ou : 

M est un vecteur de dimension (L+2) . (M+2) -4 

dont chaque composante vaut [iim. c'est a dire la valeur 
locale du module de cisaillement en chaque point discreti- 

10 se de coordonnSes (xi,ym) , ou 1 et m sont des entiers com- 
pris entre 0 et respect ivement L+1 et M+1 en eliminant les 
couples (l,m) valant (0,0), (0,M+1) , (L+1,0) et (L+1, M+1), 
L+2 et M+2 Stant les nombres de points discrStis^s dans 
1' image du milieu 1 respect ivement selon les axes X et Y, 

15 - B est un vecteur de dimension L.M, dont les 

composantes valent* -tzj^p.Uim oQ m est la pulsation de 
I'onde de cisaillement basse frequence, p est la densite 
du milieu, Ui^ est la transformee de Fourier temporelle du 
deplacement u au point de coordonnees (xi,ym)/ 1 etant 

20 compris entre 1 et L et m etant compris entre 1 et M, 

et H est une matrice de dimension L.M lignes 
sur (L+2) , (M+2) -4 colonnes, dont les composantes sont ega- 
lement toutes connues h partir de 1' equation d'onde. 

En juxtaposant un nombre suffisant d' equations 

25 (III) correspondant respect ivement a differentes disponi- 
bles du spectre de frequences de I'onde de cisaillement, 
on obtient une Equation matricielle globale qui peut etre 
resolue par inversion matricielle, pour obtenir le vecteur 
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M, c'est a dire la valeur du module de cisaillement en 
tout point du champ d' observation , 

On peut alors en deduire eventuellement la valeur 
locale de la vitesse de propagation Cs de I'onde de ci- 
saillement en chaque point, par la formule : 




Le mode de calcul serait le meme en utilisant non 
plus les deplacements, mais les deformations du milieu ob- 
serve 1 . 
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pp;yTa>jnTfTi^TTONS 

1. ProcSdfi d'imagerie utilisant les ondes de ci- 
saillement pour observer un milieu visco61ascique diffu- 
5 sant (1) qui contient des particules (5) ref lechissant les 
ondes ultrasonores de compression, procede dans lequel on 
ggngre une onde Slastique de cisaillemenc dans le milieu 
visco^lascique et on observe, au moyen d'au moins une onde 
ultrasonore de compression, le dSplacement du milieu vis- 
10 coelastique (1) soumis a ladite onde de cisaillement , 

caractSrise en ce qu'on gendre I'onde de cisaillement en 
appliquant au milieu viscoelast ique une excitation ayant 
la forme d'une impulsion basse frequence qui presente une- 
frequence centrale f comprise entre 20 et 5000 Hz, cette 
15 impulsion basse frequence presentant une duree comprise 
entre l/2f et 20/f, 

en ce qu'il comporte une etape d • observation de propaga- 
tion au cours de laquelle on observe la propagation de 
I'onde de cisaillement simultanement en une multitude de 
20 points dans le milieu observe, ces points formant un champ 
d'observation sensibleraent continu s'Stendant au moins se- 
lon un premier axe (X) , cette 6tape d'observation de pro- 
pagation de I'onde de cisaillement consistant a : 

- gmettre dans le milieu observS une succession 
25 d'au moins 10 tirs d' ondes ultrasonores de compression k 

une cadence comprise entre 100 et 100 000 tirs par se- 
conde , 

- dececter et enregistrer en temps r6el les echos 
generes par les particules ref lechissantes du milieu vis- 
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coelastique a chaque tir d'onde ultrasonore, ces echos 
correspondant a des images successives du milieu observe, 
et en ce que ledit procSdS comporte en outre une etape ul- 
t^rieure de traitement d' image au cours de laquelle on 

5 traite en temps diff^rS les images ainsi obtenues au moins 
par intercorrelation entre images successives, pour deter- 
miner en chaque point du champ d' observation un parametre 
de mouvement choisi entre le d€placement et la deformation 
du milieu viscoelastique, de faqron a obtenir ainsi une 

10 succession d' images montrant 1' evolution du parametre de 
mouvement du milieu viscoSlastique sous I'effet de la pro- 
pagation de 1 * onde de cisaillement . 

2. Procede selon la revendication 1, dans lequel 
la duree de 1' impulsion basse frequence est comprise entre 

15 l/2f et 2f. 

3. Procede selon la revendication 1 ou la revendi- 
cation 2, dans lequel la frequence centrale de 1' impulsion 
basse frequence est comprise entre 30 et 1000 Hz. 

4- Procede selon I'une quelconque des revendica- 
20 tions prec^dentes, dans lequel le milieu visco^lastique 
(1) observe est const itu6 par un corps vivant comprenant 
au moins un organe interne soumis S des mouvements impul- 
sionnels, 1' impulsion basse frequence qui gendre I'onde de 
cisaillement etant constitute par un mouvement impulsion- 
25 nel dudit organe interne. 

5. ProcedS selon l*une quelconque des revendica- 
tions 1 a 3, dans lequel le milieu viscoelastique (1) ob- 
serve est dSlimit^ par une surface exterieure (3) , et 
1' impulsion basse frequence est appliquee au niveau de 
30 cette surface exterieure. 
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6. ProcSde selon la revendication 5, dans lequel 
1' impulsion basse frequence est appliquee par un mcyen 
d' excitation choisi parmi : 

une onde acoustique genfir^e par au moins an 
5 transducteur acoustique (2) , 

et un choc gSnere localement par contact physi- 
que au niveau de la surface exterieure (3) du milieu vis- 
coelastique. 

7. Precede selon la revendication 5 ou la revendi- 
10 cation 6, dans lequel on §met les tirs d'ondes ultrasono- 

res de compression et on dStecte les Schos gSnSres par les 
particules rSf ISchissantes du milieu viscoSlast ique , au 
moyen d'une batterie de transducteurs (6) qui comprend au 
moins un transducteur (Tl, Tn) et qui est disposee au 

15 contact de la surface exterieure du milieu viscoSlast ique, 
1 ' onde de cisaillement §tant appliquSe au milieu visco- 
61astique en imposant un deplacement impulsionnel a ladite 
batterie de transducteurs . 

8. Procede selon la revendication 7, dans lequel 
20 ledit parametre de mouvement est la deformation du milieu 

viscoelast ique , 

9. Procede selon 1 ' une quelconque des revendica- 
tions precedentes, dans lequel, au cours de 1 ' observation 
de la propagation de I'onde de cisaillement, on emet entre 

25 100 et 10 000 tirs d'ondes ultrasonores de compression a 
une cadence comprise entre 1000 et 100 000 tirs par se- 
conde . 

10. Procede selon I'une quelconque des revendica- 
tions precedentes, dans lequel le champ d' observation 

30 s'etend au moins selon un plan comprenant d'une part, le 
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premier axe (X) et d* autre part, un deuxiSme axe perpendi- 
culaire au premier axe (Y) . 

11. Precede selon la revendication 10, dans le- 
quel, au cours de I'Stape d ' observation de propagation, on 

5 utilise une batterie (6) de plusieurs transducteurs acous- 
tiques (Tl, Tn) disposes au moins selon le deuxiSme axe 
(Y) pour emettre les tirs d'ondes acoustiques ultrasonores 
de compression et d^tecter les 4chos g^nSrSs par les par- 
ticules r6f ISchissantes (5) du milieu visco^lastique, les 

10 echos detectes par chaque transducteur acoustique etant 
memorises directement sans traitement prealable au cours 
de 1 • etape d ' observation de propagation, et I'etape de 
traitement d* image comprenant une sous-etape preliminaire 
de formation de voies au cours de laquelle on genere une 

15 image du milieu viscoelastique correspondant a chaque tir 
d'onde ultrasonore de compression, par combinaison d'au 
moins certains des echos re<;us par les differents trans- 
ducteurs (Tl, Tn) . 

12. Procede selon I'une quelconque des revendica- 
20 tions prScfidentes, dans lequel 1* etape de traitement 

d' images est suivie par une etape de visualisation au 
cours de laquelle on visualise au ralenti un film consti- 
tue par la succession des images traitees, chaque point de 
chaque image prSsentant un parametre opticpae qui varie 
25 suivant la valeur du parametre de mouvement affecte a ce 
point . 

13. ProcSde selon la revendication 12, dans lequel 
le paramdtre optique est choisi parmi le niveau de gris et 
le niveau chromatique. 
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14. Procede selon 1 * une quelconque des revendica- 
tions precedentes, dans lequel I'etape de traitement 
d' images est suivie par une etape de cartographie au cours 
de laquelle, a partir de 1* evolution du paramStre de tnou- 

5 vement au cours du temps dans le champ d* observation, on 
calcule au moins un parametre de propagation de I'onde de 
cisaillement en au moins certains points du champ d' obser- 
vation . 

15. Precede selon la revendication 14, dans lequel 
10 le parametre de propagation de I'onde de cisaillement qui 

est calcule au cours de I'etape de cartographie est chdisi 
parmi la vitesse des ondes de cisaillement, le module de 
cisaillement, 1 ' attenuation des ondes de cisaillement, 
l'61asticit6 de cisaillement, et la viscosite de cisaille- 
1 5 ment . 

16. Dispositif d'imagerie utilisant les ondes de 
cisaillement pour observer un milieu viscoelast ique diffu- 
sant (1) qui contient des particules (5) ref lechissant les 
ondes ultrasonores de compression, ce dispositif compre- 

20 nant des moyens d' excitation (2) pour generer une onde 
elastique de cisaillement dans le milieu viscoelastique et 
des moyens d • acquisition (CPU, Ti, Ei, Mi) pour observer, 
au moyen d ' au moins une onde ultrasonore de compression, 
le deplacement du milieu viscoelastique soumis S ladite 

25 onde de cisaillement, 

caracterise en ce que les moyens d' excitation (2) sont 
adaptSs pour appliquer au milieu viscoelastique (1) une 
excitation ayant la forme d'une impulsion basse frequence 
t qui presente une frequence centrale f comprise entre 20 et 
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5000 Hz, cette impulsion basse frequence presentant une 
duree comprise entre l/2f et 20/f, 

en ce que les moyens d' acquisition (CPU, Ti, Ei, Mi) sont 
adaptSs pour observer la propagation de I'onde de ci- 
5 saillement simultan§ment en une multitude de points dans 
le milieu observS, ces points formant un champ d' observa- 
tion sensiblement continu s'etendant au moins selon un 
premier axe (X), lesdits moyens d' acquisition 4tant adap- 
tes pour : 

10 - 6mettre dans le milieu observe une succession 

d'au moins 10 tirs d'ondes ultrasonores de compression a 
une cadence comprise entre 100 et 100 000 tirs par se- 
conde , 

dStecter et enregistrer en temps reel les echos 
15 generes par les particules ref lechissantes du milieu vis- 
coelastique k chaque tir d'onde ultrasonore, ces echos 
correspondant a des images successives du milieu observe, 
et en ce que ledit dispositif comporte en outre des moyens 
de traitement d» image {CPU, S, DSP) adaptes pour traiter 
20 en temps differe les images obtenues 1 partir des moyens 
d • observation , au moins par intercorrSlation entre images 
successives, pour determiner en chaque point du champ 
d' observation un paramStre de mouvement choisi entre le 
dSplacement et la deformat ion du milieu viscoelastique , de 
25 fagon S obtenir ainsi une succession d' images montrant 
1' evolution du parametre de mouvement du milieu viscoelas- 
tique sous I'effet de la propagation de I'onde de ci- 
saillement . 

17. Dispositif selon la revendication 16, dans le- 
30 quel les moyens d' observation comprennent une batterie 16) 
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de transducteurs (Tl, Tn) qui inclut au moins un trans- 
ducteur et qui est adapt^e pour etre disposSe au contact 
d'une surface exterieure (3) dSlimitant le milieu visco- 
elastique (1), les tnoyens d'excitation (2) etant adaptes 
pour imposer un deplacement impulsionnel a ladite batterie 
(6) de transducteurs. 

18. Dispositif selon la revendication 11, dans le- 
C[uel ledit parametre de mouvement est la deformation du 
milieu viscoelastique . 

19. Dispositif selon 1 ' une quelconque des revendi- 
cations 17 et 18, dans lequel le champ d' observation 
s'etend au moins selon un plan comprenant d'une part, le 
premier axe (X) et d* autre part, un deuxieme axe (Y) per- 
pendiculaire au premier axe, la batterie ^ (6) de transduc- 
teurs comportant plusieurs transducteurs (Tl, Tn) dis- 
poses au moins selon le deuxieme axe, des moyens de com- 
mande (4, CPU) Stant prfivus pour s61ectivement faire fonc- 
tionner le dispositif soit dans le mode d'imagerie par on- 
des de cisaillement , soit dans un mode d ' echographie stan- 
dard permettant d'acqu^rir entre 10 et 100 images par se- 
conde . 
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